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ABSTRACT 

Some of the satellites in the Solar System, including the Moon, appear to have been captured from heliocentric 
orbits at some point in their past, and then have evolved to the present configurations. The exact process of how 
this trapping occurred is unknown, but the dissociation of a planetesimal binary in the gravitational field of the 
planet, gas drag, or a massive collision seem to be the best candidates. However, all these mechanisms leave the 
satellites in elliptical orbits that need to be damped to the present almost circular ones. Here we give a complete 
description of the secular tidal evolution of a satellite just after entering a bounding state with the planet. In 
particular, we take into account the spin evolution of the satellite, which has often been assumed synchronous 
in previous studies. We apply our model to Triton and successfully explain some geophysical properties of this 
satellite, as well as the main dynamical features observed for the Neptunian system. 

Subject headings: celestial mechanics — methods: analytical — planets and satellites: general — planets and 
satellites: individual (Neptune, Triton) 



1. INTRODUCTION 
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■ Both the Earth's Moon and Pluto's moon, Charon, have an 
'"'important fraction of the mass of their systems, and therefore 
Q they could be classified as double-planets rather than as satel- 
$^ lites. The proto-planetary disk is unlikely to produce such 
Systems, and their origin seems to be due to a catastrophic 
Ci jmpact of the initial planet with a body of comparable dimen- 
'""'sions (e.g. Canup and Asphaug 2001; Canup 2005). On the 
^ bther hand, Neptune's moon, Triton, and the Martian moon, 
t'hobos, are spiraling down into the planet, clearly indicating 
JlC that the present orbits are not primordial, and may have un- 
I dergone a long evolving process from a previous capture (e.g. 
^ Mignard 1981; Goldreich et al. 1989). 

i The present orbits of all these satellites are almost circular, 
and their spins appear to be synchronous with the orbital mean 
motion, as well as being locked in Cassini states (e.g. Colombo 
1966; Peale 1969). This also applies to the Galilean satellites 
of Jupiter, which are likely to have originated from Jupiter's 
accretion disk and additionally show orbital mean motion res- 
onances (e.g. Yoder 1979). All these features seem to be due 
to tidal evolution, which arises from differential and inelastic 
deformation of the planet by a perturbing body. 

Previous long-term studies on the orbital evolution of satel- 
lites have assumed that their rotation is synchronously locked, 
and therefore limits the tidal evolution to the orbits (e.g. Mc- 
Cord 1966). However, these two kinds of evolution cannot be 
dissociated because the total angular momentum must be con- 
served. Additionally, it has been assumed that the spin axis 
is locked in a Cassini state with very low obliquity. Although 
these assumptions are correct for the presently known situa- 
tions, they were not necessarily true throughout the evolution. 

The aim of this Letter is to model the orbital evolution of 
a satellite from its origin or capture until the preset day, in- 
cluding spin evolution for both planet and satellite, and also 
to make predictions regarding its future evolution. We provide 
a simple averaged model adapted for fast computational sim- 
ulations, as required for long-term studies. In the last section 
we apply this model to the Triton-Neptune system. 
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2. THE MODEL 

We consider a hierarchical system composed of a star, a 
planet and a satellite, with masses M ^ mp ^ mi, respec- 
tively. Both planet and satellite are considered oblate ellipsoids 
with gravity field coefficients given by and J21 , rotating 
about the axis of maximal inertia along the directions Sq and 
§1, with rotation rates ujq and wi, respectively. The potential 
energy U of the system is then given by (e.g. Smart 1953): 
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where terms in [Ri/rjY have been neglected (i, j = 0, 1). G is 
the gravitational constant, Ri the radius of the planet or the 
satellite, r; the distance between the planet and the star or the 
satellite, and P2{x) = {ix^ — l)/2 the Legendre polynomial of 
degree two. 

Neglecting tidal interactions with the star, the tidal potential 
is written (e.g. Kaula 1964): 
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where fc2i is the potential Love number for the planet or the 
satellite, and t[ the position of the interacting body at a time 
delayed of Atj. For simplicity, we will adopt a model with 
constant Aii, which can be made linear (e.g. Mignard 1979; 
Hut 1981; Neron de Surgy and Laskar 1997): 



r'; ~ ri -f Mi (wiSj X ri - fi) 



(3) 



The complete evolution of the system can be tracked by the 
evolution of the rotational angular momentums, ~ CiLUiSi, 
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the orbital angular momentums, ~ miniaf{l — e^)^/^ki, 
and the orbital energy Ei ~ —Gmomi/{2ai). rii is the mean 
motion, aj the semi-major axis, the eccentricity, and Cj the 
principal moment of inertia. The contributions to the orbits 
are easily computed from the above potentials as 

Lo = -roxVr„Z7; Li = -nxVr.C/; Ei = -h-^r^U . (4) 

Since the total angular momentum is conserved, the contri- 
butions to the spin of the planet and satellite can easily be 
computed from the orbital contributions: 



Ho + Hi + Lo + Li = . 



(5) 



Because tidal effects act in long-term time-scales, wc further 
average the equations of motion over fast angles, namely the 
true anomaly and the longitude of the periapse. The resulting 
equations for the conservative motion are (e.g. Boue and Laskar 
2006): 
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and 



cos£i = s, • ko , cos^i = Sj • ki , cos / = ko • ki , (11) 

are the direction cosines of the spins and orbits: Si is the obliq- 
uity to the orbital plane of the planet, 6i the obliquity to the 
orbital plane of the satellite, and I the inclination between 

orbital planes. 

For the dissipative tidal effects, we can obtain the equations 
of motion directly from equation (4), using Ut instead of U, 

that is, Lo = 0, 

TT ( f ( xS*+COS^ikl - \ 

Hi = -Kim /i(ei) /2(ei)ki , (12) 

\ 2, m I 

Ei=Y,Kinl(f2{ei)cos0i—-Mei)) , (13) 



where = 0, 1), 

Ki = Ati{3k2,Gm''jR^i)/at , (14) 

and /i(e) = (l + 3e2 + 3eV8)/(l-e2)9/2 , /2(e) = (1-f 15eV2-f 
456^^/8 + 5eVl6)/(l - e^)^ , /3(e) = {1 + 31eV2 -|- 255eV8 + 
185eVl6 + 25eV64)/(l - e^^^/^ . 



3. SECULAR EVOLUTION 

In the previous section we presented the equations that rule 
the tidal evolution of a satellite in terms of angular momenta 
and orbital energy. However, the spin and orbital quanti- 
ties are better represented by the rotation angles and ellip- 
tical elements. The direction cosines (Eq.ll) are obtained 
from the angular momenta vectors, since Sj = Hj/||Hj|| and 
ki = Li/||Li||, as well as the rotation rate cui = Hi -Si/Ci. The 
semi-major axis and the eccentricity can be obtained from Ei 
and ||Li||, respectively. 

The variation in the satellite's rotation rate can be computed 
from equation (12) as = Hj • Sj/Cj, giving (Correia and 
Laskar 2009): 
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For a given obliquity and eccentricity, the equilibrium rotation 
rate, obtained when tbi = 0, is attained for: 

t^i ^ /2(ei) 2cos6ii 

ni /i(ei) 1 + 0082^1 ' ^ ^ 

The obliquity variations can be obtained from equation (11): 
dcosOj _ ilj ■ (ki - cosgjSj) ^ Li ■ (§» - cosgjki) ^^^^ 
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For the conservative motion (Eqs.6, 7), stable configurations 
for the spin can be found whenever the vectors (§1, ki, ko) 
or (§1, ki, §0) are coplanar and precess at the same rate g 
(e.g. Colombo 1966; Peale 1969). The first situation occurs if 
7 ^ /3o (outer satellite) and the second situation when 7 <^ /3o 
(inner satellite) . The equilibrium obliquities can be found from 
a single relationship (e.g. Ward and Hamilton 2004): 
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where Ai = Pi/{CiU)ig) is a dimensionless parameter and Iq 
is the inclination of the orbit of the satellite with respect to 
the Laplacian plane {Iq — I and g ~ 7cos//||Li|| for an outer 
satellite, and Iq — Oo and g c± /3o cos^o/||iii|| for an inner 
satellite) (e.g. Laplace 1799; Mignard 1981; Tremaine et al. 
2009). The above equation has two or four real roots for Oi, 
which are known by Cassini states. In general, for satellites 
we have Jq ~ 0, and these solutions are approximately given 
by: 
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For a generic value of Jo, when Ai ^ 1, which is often the 
case of an outer satellite, the first expression gives the only 
two real roots of equation (18), one for 61 ~ Iq and another for 
9 1 ~ tt—Iq- On the other hand, when Ai ^ 1, which is the case 
of inner satellites, we will have four real roots approximately 
given by expressions (19). 

In turn, the dissipative obliquity variations are computed by 
substituting equation (12) in (17) with ||Hi|| -C ||Li||, giving: 
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Because of the factor ni/wi in the magnitude of the obliquity 
variations, for an initial fast rotating satellite, the time-scale 
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for the obliquity evolution will be longer than the time-scale for 
the rotation rate evolution (Eq.l5). As a consequence, it is to 
be expected that the rotation rate reaches its equilibrium value 
(Eq.l6) earlier than the obliquity. Thus, replacing equation 
(16) in (20), we have: 
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We then conclude that the obliquity can only decrease by tidal 
effect, since 0i < 0, and the final obliquity tends to be captured 
in low obliquity Cassini states. 

The variations in the semi-major axis are obtained from the 
energy variations Ei , 



^1 = X] i -^2(^1) cos6'j^ - /3(ei) 
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while the eccentricity is obtained from the norm of the orbital 
angular momentum ||Li|| = miniof (1 — ef)^/^: 

miaj \ 18' 



-2 \ 7^/4(61) cos6'i^ - /5(ei) ) ei 



(23) 



where /4(e) = (1 + 3eV2 + eV8)/(l - e^)^ , /5(e) = (1 + 
15eV4 + 15e4/8 + 5eV64)/(l - e^^/^ . For gaseous planets 
and rocky satellites we usually have k2gAto ^ Ati, and we 
can retain only terms in Ki . 

The ratio between orbital and spin evolution time-scales 
is roughly given by Ci/(miaf) <C 1, meaning that the spin 

achieves an equilibrium position (Hi = 0) much faster than 
the orbit. Replacing the equilibrium rotation rate (Eq. 16) with 
^1=0 (for simplicity) in equations (22) and (23), gives: 

ai = -^/6(ei)e? , (24) 
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(25) 



where /6(e) = (1 + 45eVl4 + 8e* + 685eV224 + 2556^448 + 
256^71792) (1 - e2)-i5/2/(i + 3^2 _^ '^pi/gj . Thus, we always 
have ai < and ei < 0, and the final eccentricity is zero. 
Another consequence is that Li — — Hq ~ 0, and the quantity 
ai(l — ef) is conserved. The final equilibrium semi-major axis 
is then given by a/ = ai(l — ef) . However, from this point 
onwards, the tidal effects on the planet cannot be neglected 
(Eq.22), and they govern the future evolution of the satellite's 
orbit. For af < as or 0q > 7r/2, where = Gmo / {ujq cos 0q)'^ , 
the semi-major axis continues to decrease until the satellite 
crashes into the planet, while in the remaining situations it 
will increase. 

4. APPLICATION TO TRITON-NEPTUNE 

Neptune's main satellite, Triton, presents unique features 
in the Solar System. It is the only moon-sized body in a 
retrograde inclined orbit and the images taken by the Voy- 
ager 2 spacecraft in 1989 revealed an extremely young surface 
with few impact craters (e.g. Cruikshank 1995). This satel- 
lite should have remained molten until about 1 Gyr ago and 
its interior is still warm and geologically active considering its 
distance from the Sun (Schenk and Zahnle 2007). Its com- 
position also presents some similarities with Pluto (Tsurutani 
et al. 1990). 



These bizarre characteristics lead one to believe that Triton 
originally orbited the Sim, belonging to the family of Kuiper- 
belt objects. Most likely during the outward migration of Nep- 
tune, the orbits of the two bodies intercepted and capture oc- 
curred. This possibility is strongly supported by the fact that 
Triton's present orbit lies between a group of small inner pro- 
grade satellites and a number of exterior irregular satellites 
both prograde and retrograde. Nereid, with an orbital eccen- 
tricity around 0.75, is also believed to have been scattered from 
a regular satellite orbit (McKinnon 1984). 

How exactly the capture occurred is still unknown, but some 
mechanisms have been proposed: gas drag (Pollack et al. 1979; 
McKinnon and Leith 1995), a collision with a pre-existing reg- 
ular satellite of Neptune (Goldreich et al. 1989), or three-body 
interactions (Agnor and Hamilton 2006; Vokrouhlicky et al. 
2008). All these scenarios require a very close passage to Nep- 
tune, and leave the planet in eccentric orbits that must be 
damped by tides to the present one. Tides are thus the only 
consensual mechanism acting on Triton's orbit. The tidal dis- 
tortion of Triton after a few close passages around Neptune, 
and the consequent dissipation of tidal energy, can account 
for a substantial reduction in the semi-major axis of its orbit, 
quickly bringing the planet from an orbit outside Neptune's 
Hill sphere (~ 4700 i?o) to a boimded orbit. Therefore, it can- 
not be ruled out that Triton was simply captured by tidal 
interactions with Neptune after a close encounter in an almost 
parabolic orbit (McCord 1966). 

In this Letter we simulate the tidal evolution of the Triton- 
Neptune system using the complete model described in Sect. 2. 
Triton is started in a very elliptical orbit with ei = 0.9968 
and a semi-major axis of oi = 2354 i?o, corresponding to a 
final equilibrium ay ~ 15i?0) close to the present position of 
14.33 i?o. These specific values also place the satellite outside 
the Hill sphere of Neptune at the apoapse and give a closest 
approach at a periapse of 7.53 i?o- However, the non-secular 
perturbations of the Sim on Triton's orbit will cause the eccen- 
tricity to vary around the mean value, allowing the periapse 
and apoapse distances to attain lower and higher values (Gol- 
dreich et al. 1989). 

For the radius of the bodies we use Rq = 24 764 km and 
i?i = 1353km (Thomas 2000), while for the masses, the J2 
of Neptune and the remaining orbital and spin parameters we 
take the present values as determined by Jacobson (2009). For 
Triton we adopt J2 = 4.38 x 10""', the value measured for Eu- 
ropa (Anderson et al. 1998), and C22 = 0, since our model 
does not take into account spin-orbit resonances. This choice 
is justified because Triton's observed topography never varies 
beyond a kilometer (Thomas 2000). In addition, Triton shoidd 
have undergone frequent collisions either with other satellites 
of Neptune, or with external Kuiper-belt objects, and any cap- 
ture in a spin-orbit resonance different from the synchronous 
one, may not last for a long time (Stern and McKinnon 2000). 
For tidal dissipation we adopt the same parameters as previ- 
ous studies, that is, = 0.407 and Qo = 9000 (Zhang and 
Hamilton 2008), and fca^ = 0.1 and Qi = 100 (Goldreich et al. 
1989; Chyba et al. 1989), where Qr^ = coiAU. 

As for the orbit, the initial spin of Triton is unknown. We 
tested several possibilities, but tides acting on the spin always 
drive it in the same way: the rotation rate quickly evolves into 
the equilibrium value given by equation (16), while the obliq- 
uity is trapped in Cassini state 2, for 61 ~ 180° — /q = 23° 
(Eq.l9). In our standard simulation (Fig.l) we start Triton 
with a rotation period of 24 h and an obliquity 9i = 170° (li- 
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Fig. 1. — Secular tidal evolution of the Triton-Neptune system. We 
plot the orbital inclination to the Laplacian plane Iq, the obliquity of 
Neptune to the ecliptic £q and the obliquity of Triton to the orbital plane 
6\ (top), and the semi- major axis ratios ai/Ro and af/Ro, the inverse of 
the eccentricity e^"*^ and Triton's rotation rate ratio oJi/ni (bottom). 

brating around Cassini state 3). In the very beginning, only 
Cassini states 2 and 3 exist (Ai <C 1), but according to equa- 
tion (21), equilibrium in Cassini state 3 cannot be maintained 
(Fig. la). As the semi- major axis decreases, the equilibrium 
obliquity for state 2 increases toward 90° (Eq.l9). At some 
point, the tidal torque becomes stronger than the conserva- 
tive torque and the spin quits this state. The obliquity subse- 
quently evolves into Cassini state 1 with 6i ~ Iq/Xi (Eq.19). 
The orbital inclination to the Laplacian plane is more or less 



constant, but it presents some reduction when tides raised 
by Triton on Neptune become important (present day situ- 
ation). At the very end of the evolution (ai ~ 7.79 i?o), this 
trend is reversed and the inclination quickly evolves into 180°. 
The obliquity of Neptune (£o) does not undergo any signif- 
icant dissipation, but presents a secular oscillation of about 
one degree, from the moment Triton becomes an inner satel- 
lite (ai < 100 i?o)- The semi-major axis and the eccentricity 
always decrease, as predicted by equations (24) and (25), and 
the quantity a/ = ai(l — e^) is preserved, during the first stages 
of the evolution, with the reduction observed being caused by 
tides on Neptune. The eccentricity is very high during the first 
evolutionary stages, decreasing very slowly. However, as soon 
as Triton becomes an inner satellite, the eccentricity quickly 
drops to a value very close to zero. Finally, the rotation of the 
satellite decreases as the satellite orbit shrinks into Neptune. 
It presents a rapid decrease when Cassini state 2 approaches 
90°, explained by Eq. (16), and ultimately stabilizes in the 
synchronous resonance, the presently observed configuration. 
5. CONCLUSION 

We have shown that the tidal evolution of a satellite can 
only be correctly modeled when its spin is taken into account. 
Previous studies adopted synchronous motion from the very 
beginning, which is not the case for eccentric orbits. Tidal dis- 
sipation was therefore underestimated. With the same tidal 
parameters used by Goldreich et al. (1989) we are able to cir- 
cularize Triton's orbit in only 1 Gyr. Different tidal models 
and parameters may change the time-scale for the evolution, 
but the global picture should remain the same. We are also 
able to explain the present small value of Triton's eccentricity 
ei ~ 10~^, as well as its small obliquity 9i = 0.46° and past 
evolution through Cassini states. In particular, we can exclude 
the possibility that Triton was initially captured in state 1 and 
determine when exactly the transition from state 2 to state 1 
occurred (Chyba et al. 1989). 

Our study should also apply to the Moon, Charon and the 
satellites of Mars, although in this case we need to take into 
account the quadropole moment of inertia C22 7^ (Correia 
2006). It can also be easily generalized to other stellar systems. 

The author thanks J. Laskar for helpful discussions. This work 
is dedicated to the memory of Antonio and Nazareth Morgado. 
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